Mean-Field Approach to Charge, Orbital, and Spin Ordering in Manganites 
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We present a mean-field theory of charge, orbital, and spin ordering in manganites at 50% and 
0% dopings by considering Jahn- Teller interaction, nearest-neighbor repulsion, and no single-site 
double occupancy. For spinless fermions, we show that Jahn- Teller distortion and charge-orbital 
ordering occur simultaneously. In our two-dimensional model at 50% doping, for small nearest- 
neighbor repulsion the system is orbitally polarized while for larger repulsion the system undergoes 
CE type ordering. As for the 0% doping case, the ground state is orbitally antiferromagnetic. Upon 
including spin degree of freedom, at both 50% and 0% dopings the ordering remains the same at 
small antiferromagnetic coupling between adjacent core spins. 
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Manganites of the form A\_ x B x MnOz (with A=La, 
Pr, etc; B=Ca, Sr, etc) display a rich variety of phenom- 
ena and phases as the doping x is varied Q. In particu- 
lar close to 0% doping the systems display orbital order- 
ing at higher temperatures while at lower temperatures 
magnetic ordering results 0. In narrow band materi- 
als like Pra.sCao.sMnOa, as the temperature is lowered 
to Tco the system undergoes CE type charge-orbital 
(CO) ordering and at a much lower temperature TV pe- 
culiar antiferromagnetic (AF) order sets in (see Fig. 
[I]). In wider band materials like Nd , 5 Sr , 5 MnO 3 fer- 
romagnetic ordering occurs at higher temperatures and 
at lower temperatures simultaneous charge, orbital, and 
spin (AF) ordering results ||[|. In layered-perovskites 
like Lao^Sri.sAInOi it has been reported that Tco is 
higher than Tjv In systems where Tco and Tjv 

coincide, Solovyev and Terakura argue that the CO- 
ordering is a result of the peculiar AF ordering along 
with Hund's rule physics 0. van den Brink et al. too 
analyze systems like Nd a 5 Sr ^MnO^ and find that be- 
cause of a topological factor in the hopping, the one- 
dimensional ferromagnetic zigzag chains are orbitally or- 
dered However Refs. $M do not include electron- 
phonon interaction physics nor do they consider narrow 
band systems like Pr ^CaQ.^MnO 3 . Evidence of sig- 
nificant change in lattice parameters has been reported 
at Tco for Pn- x Ca x Mn0 3 (at x=0.4, 0.5) @ and in 
La . 5 Ca . 5 MnO 3 @. 

In this paper we derive a general framework for an- 
alyzing the concomitant charge density wave (CDW) 
and orbital density wave (ODW) instability by consid- 
ering Jahn- Teller interactions and showing that the two- 
dimensional (2D) ordering can be understood in terms of 
the one-dimensional (ID) Kohn anomaly at wavevector 
q = 2k p. The problem essentially reduces to comparing 
the energy resulting from the hopping term for various 
possible ordered states. We show that Jahn- Teller dis- 
tortion is accompanied by CO-ordering for the carriers. 
We find that above a critical value of nearest-neighbor 



repulsion and for small antiferromagnetic coupling be- 
tween adjacent core spins, CE type charge, orbital, and 
spin (COS) ordering is realized at 50% doping. At % 
doping orbitally antiferromagnetic order results. 

We first study the ordering phenomena without spin 
degree of freedom because Tco > Tn- We begin by con- 
sidering CDW instability for a ID Holstein model given 
below: 



H = 



(ij) <f 



(i) 



q, k 



where a,j (fg-) is the electron (phonon) destruction op- 
erator, (ij) corresponds to nearest neighbors, loq is the 
optical phonon frequency (h = 1), A(< 0) is the electron- 
phonon coupling, t is the hopping integral, and N is the 
number of sites. We know that in ID the non-interacting 
polarizability xo (o) (Lindhard polarizability) diverges at 
wavevector q = 2kp. Then even in the weak-coupling 
limit we will have a CDW instability because the ex- 
pectation value (X)fc a | + - a fe) diverges in the interaction 
part of the Hamiltonian. Furthermore, upon calculat- 
ing the double derivative with respect to time of the 
ionic position coordinate Qg- = y/l/(2MuJo){fq + f-$) 
by using the relation = —[[Qq,H],H] we get Qq = 
—^oQq ~ 2^4 2 woXo( ( z) ( 9g- Thus we see that the renor- 
malized phonon frequency becomes soft even for vanish- 
ingly small electron-phonon interaction leading to lattice 
distortion. The order parameter for the CDW state is 
(fj) = Aje l ^ j . Then within mean-field we get from Eq. 
(|), after averaging over phonon coordinates, the follow- 
ing (see Ref. fllj for CDW in metals): 

H = -t ^2 a \ a i + A| + 2A^2 a ] a j^j cos <t>j ■ ( 2 ) 

(ij) i 3 

The energy in the above equation is minimized for cf>j = 
and for u>oAj = — Afajdj). 
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We will now solve for the CDW state at half-filling, 
using Aj = A /2 + Ai cos(7rj). We first note that A = 
—A/ujq and that Ai < Ao/2. Next we add the nearest- 
neighbor interaction term Vnn = V J2j s n j n j+s with 
S being nearest-neighbors and make the mean-field ap- 
proximation rijUj+i = nj(rij + i) +nj + i(rij) — (n,j+i){nj). 
Then the Hamiltonian in momentum space, on folding 
the Brillouin zone, gets modified to be 

H = Y} Hk + Hk +- + 2AB+A 1 (aj[a fc+7r + H.c.)] 
k 

+uj N{B^A 2 /4 + B+A 2 1 ), (3) 

where Hk = (— 2icos/c + AB- Ao)atafc, B± = 1 ± 
zVujo/A 2 with z being the coordination number, and the 
reduced Brillouin zone is — 7r/2 < k < tt/2. Then diago- 
nalizing the above Hamiltonian yields the eigen energies 

E± = AB_A ± ^4t 2 cos 2 k + AA 2 B\A\. Upon mini- 
mizing the total energy with respect to Ai we get the 
optimum A x » -A/(2w ) [l - (A/B+) 2 + 1/4(A/B+) 4 ] 
for small values of X/ B + where the polaron size parame- 
ter A = tujo/A 2 . Then the minimized energy per particle 
is -A 2 /lu [l + X 2 /B + - A 4 /(4S|)] where the last two 
terms involving A represent the hopping term energy. We 
see that including nearest-neighbor interactions does not 
alter the results qualitatively although it enhances the 
charge modulation. 

We will now proceed to consider manganite systems 
with two orbitals per site. The Hamiltonian consists 
of three parts H = Hi + H 2 + H3 where Hi is the 
hopping term, H2 the ionic term, and H$ the electron- 
ion interaction term. The hopping term is given by 
—t y)uj\ e l e j where (ij) corresponds to nearest-neighbors 
in the a(— x,y, or z) direction and e is the destruction 
operator of the orbital ip 3a 2__ r 2. Here it should be men- 
tioned that for propagation in the x-direction the orbital 
basis is ip^ x 2_ r i and ip y 2_ z 2, for y-direction it is ^3y 2 -r 2 
and ipz 2 -x 2 j etc. 

To analyze the ordering problem, we consider 
the orthonormal wavefunctions ipx[Y] — ^gKv** + 

l)V'3x 2 [j/ 2 ]-r' 2 + (v3 — l)V'3a 2 b; 2 ]-r 2 ] the justification for 
which will be given below |12| . The interaction term in 
this basis is given by 

H 3 = A^2Mv y^iQ2j(a X i a xj - a Yj a Y j) 
3 

+Q3j( a xj a Yj + a Yj a X j)], (4) 

where a X j and a Y j are the destruction operators of the 
electrons in ipx and ipY orbitals at site j, with Q2 and 
Q3 corresponding to the two normal modes. The ionic 
part of the Hamiltonian is then given by 

^ = f + <4o + f + Q%)- ( 5 ) 
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We further enforce that at most only one orbital can be 
occupied at each site by setting {a X j a xj) ( a Yj a Yj) = 0. 

We study the problem in 2D (say the xy-plane) and 
observe that the overlap between ipx an d i/'3x 2 -r 2 is 
greater than 96.6% and (ipx^p3x 2 -r 2 } — (^/VVV-z 2 ) = 
(4>Yip3 y 2 -r 2 } — \(4 l xipz 2 -x 2 )\- Hence based on the hop- 
ping term of the Hamiltonian, to bring out the essential 
physics, we approximate Hi to be 

Hi re -2t ( cos k x a xk a xk + cos k v a Yj : a Yk) ■ ( 6 ) 
fe 

The 50% doping case will be first considered. In H3 
we note that (VVa^ - a v r — a} -. a,,r) diverges at 

\L^h Xk+q Xk Yk+q «' ° 

q = Q/2 = [tt/2, tt/2] for the orbitally unpolarized case 
and at q = [it,...] or q = [...,7r] for the orbitally po- 
larized case leading to a concomitant CDW and ODW 
instability with the phonon frequency corresponding to 
Q2q going soft. To see this transparently we first recog- 
nize that for the electrons in the ipx (ipy) orbitals the 
dispersion relation is like the ID case due to which the 
Fermi sea is rectangular with — ir/4 < k x (k y ) < tt/A 
and —it < k y (k x ) < it for the orbitally unpolarized 
case. Hence, the Lindhard polarizability diverges at 
[2k X Fi ■■•] and [...,2k y F] for ipx and ipy electrons respec- 
tively. As for the orbitally polarized case, on choos- 
ing ipy orbitals without loss of generality, the Fermi 
sea is given by —tt/2 < k y < tt/2; —it < k x < tt 
and thus the non-interacting polarizability is singular 
at [2k yF ,...\. Furthermore %+ Ji x% - a f y - +f a rf ) 

((Efc a X k +? a Yk + a Y k + q a xk^ is similar to the longitudi- 
nal (transverse) spin susceptibility in a Hubbard model 
with ipx and ipy orbitals standing for up and down spins. 
Thus the system, to lower its energy, prefers H3 to be di- 
agonal in the ipx and ipy orbital basis with only Q2 mode 
getting excited at all sites. The Hamiltonian after aver- 
aging over the phonon coordinates yields the following 
two-orbital 2D version of Eq. (^|) with optimum values 
of <pj and Aj\ 

H = -2tJ2 ( cos k ^ a xi a xk + cos k v a Y % a Y%) 

k 

A 2 

-2— Yl( a Xj a xj - a Yj aY ){a X] a Xj - a Yj a Y j) 
3 

A 2 

+ — Y.{ a xj a Xj ~ a Yj a Y j) 2 . (7) 

3 

We will now compare the energies for the CDW states 
in the orbitally polarized and unpolarized states. For the 
orbitally unpolarized case the order parameter (a>xj a Xj ~ 
a Yj a Yj) = Ccos(<5/2 • Rj). Then using the fact that 

( a Xj a Xj)( a Yj a Yj) — 0, we see in the present 50% dop- 
ing case that \C\ = 1 and that we get CE type CO- 
ordering. Then the CDW order parameter is given by 
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(a X j a Xj + ay-ayj) — cos 2 (Q/2 • Rj). We next intro- 
duce nearest- neighbor interaction Vnn arL d note that 
(nj){rij+s) = and thus the CDW is also compatible 
with nearest-neighbor repulsion. Then taking C = — 1 
without loss of generality, the total Hamiltonian is 



,t 



H = 



k,a=X,Y 



£ ^k+nQ/2 + 
n=0 



A 2 N 

2oj n 



(8) 



is over the reduced Brillouin zone — 7r/4 < k x ( y \ < 7r/4; 



-7r < < 7r and 



We will now include spin degree of freedom and analyze 
the magnetic ordering that results from the CO-ordering 
derived above. We consider antiferromagnetic spin cou- 
pling between adjacent localized core spins of the form 
i s Sj ' Sj+S (J > and S = 3/2) and an infinitely 
strong Hund's coupling between itinerant electrons and 
localized spins. Then in the Hamiltonian, the kinetic en- 
ergy becomes -tJ2uj] cos(6ij /2)(a Xi axj+Oy i aYj) where 
[ij] are nearest-neighbors in the x(y)-direction for ipx(Y) 
where the momentum summation, for ipx(Y) electrons, orbital electrons, % is the angle between adjacent core 

spins Si that are treated classically. As for the 2D two 
orbital case at 50% doping, for a CE type magnetic or- 
dering we get the same expectation value for the kinetic 
energy because 0y = when an electron in the ipxii^y) 
orbital tries to jump to the next unoccupied orbital in the 
x- (y-) direction (see Fig. |l|). The nearest-neighbor core 
spin coupling gives zero value. The contribution from the 
remaining spin independent terms in the Hamiltonian is 
the same as before and thus the total energy per particle 
remains the same for this spin order. On the other hand 
if the system is totally spin antiferromagnetic, the kinetic 
energy is zero and the core spin interaction energy per 
particle is — AJS 2 . Thus we see that if the magnitude of 
the hopping term energy per particle is more than 4 JS 2 
we get CE type COS ordering in 2D p). As for the 



H X(Y) 



[-2tcosk x + Vz)a} x{Y)% a x{Y)% 

-( + ) AA ^ a x(Y)k a X(Y)k+Q/2 + H - C - 

-Vz/2\a) , - *a Yfv sj> + H.c.]. 

1 1 X(Y)k+Q X(Y)k J 



(9) 



We diagonalize the Hamiltonian in Eq. (||) to obtain 

Vx(Y) + Wzrfx {Y) + 4 V 2 x{Y) [(V 2 z 2 - A 2 A 2 ) - t 2 ] 
-8rix(Y)Vzt 2 + I6t 2 sm 2 (k x(y) )cos 2 (k x(y) ) = 0, (10) 



where rjx{Y) = 



X(Y) 



2Vz with E- 

k 



X(Y) 



being the eigen 



energies for i/jx(y) electrons. Treating A as a small pa- 
rameter we obtain after some algebra the following ex- 
pressions for the lowest energies 



E* {Y) = 



u> 



X 4 



s^ 2 (2k x{y) ) + — - 



where terms of order t\ 5 or higher are neglected. On fill- 
ing up the lowest band we obtain the energy per particle 
to be -A 2 /u Q [l + X 2 /B+ - X 4 /{ABl )(l/2 + B-/B+)]. 

As for the orbitally polarized case, the order param- 
eter is given by Aj = A /2 + Ai cos(Q • Rj) and the 
problem is similar to the ID one-orbital case Jl^|. On 
choosing the ipy orbitals to be occupied, the reduced 
Brillouin zone is —tt/2 < k y < tt/2; —tt < k x < tt. 
Then the minimized energy per particle is, as before, 
-A 2 /lu [1 + X 2 /B + - X A /(AB\ )]. It is of interest to note 
that for small nearest-neighbor repulsion V, the energy 
of the orbitally polarized state is lower than that of the 
orbitally unpolarized state while for Vz > A 2 /uiq it is 
higher than that of the CE type state. 

We will now consider the 0% doping case. Here 
we get an ODW instability at q = Q leading to an 
orbitally antiferromagnetic state. The reduced Bril- 
louin zone for i^xCY) orbital is — 7r/2 < k x {k y ) < 7r/2; 
— 7r < k y (k x ) < tt. The ground state energy per par- 
ticle when A is a small parameter is obtained to be 
Vz - A 2 /lu q [1 + X 2 /2 - 3A 4 /32]. As for the orbitally 
polarized case, the energy per particle is — A 2 /ujq + Vz. 
Thus we see that the energy for the orbitally unpolarized 
case is lower than that for the orbitally polarized case at 
0% doping. 



case when the system is orbitally polarized (say of tpy 
orbitals) , if the magnetic ordering comprises of spin po- 
larized chains in the y-direction that are antiferromagnet- 
ically coupled in the x-direction (see Fig. ||), the energy 
per particle is still the same as when the spin effects are 
ignored [^5| . On the other hand if the orbitally polarized 
case is completely spin antiferromagnetic again the ki- 
netic energy is zero and the core spin interaction energy 
per particle is — AJS 2 . Thus here too we see that if the 
magnitude of the hopping term energy per particle is less 
than 4JS 2 , AF order results (see Fig. ||). 

As for the 0% doping case, when the magnitude of 
the hopping term energy per particle is more than AJS 2 , 
similar analysis as above yields an orbitally antiferromag- 
netic and spin ferromagnetic state. Otherwise, the sys- 
tem is spin antiferromagnetic. 

The above analysis can be extended to other dopings as 
well. All in all, we feel that our analysis explains the COS 
ordering in layered compounds like Lao.5SV1.5A/nO4 
and narrow-band systems like Pro.5Cao.5Af 71O3 where 
Tco > Tn. Since we have shown that magnetic ordering 
results from the CO ordering, our analysis should also 
hold for wider-band systems like NdQ^Sr^.^MnO^. 
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FIG. 1. CE ordering in x-y plane. Exaggerated Q2 distortion is shown by the dash-dot line and the zig-zag ferromagnetic 
chains by dashed lines. The Mn 3+ sites are depicted by $x,y orbitals and core spins, the Mn 4+ sites by only core spins, and 
the oxygen sites by O- 

FIG. 2. Orbitally polarized state in x-y plane. Exaggerated Q2 distortion is shown by the dash-dot line and the ferromagnetic 
chains by dashed lines. The Mn i+ sites are depicted by $y orbitals and core spins, the Mn 4+ sites by only core spins, and 
the oxygen sites by O- 



FIG. 3. Phase diagram for the 50% doping case in 2D. 
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